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Fractional photon-assisted tunnelling is investigated both numerically and analytically in a double- 
well lattice. While integer photon- assisted tunnelling is a single-particle effect, fractional photon- 
assisted tunnelling is an interaction- induced many-body effect. Double- well lattices with few parti- 
cles in each double well are ideal to study this effect far from the mean- field effects. It is predicted 
that the 1/4-resonance is observable in such systems. Fractional photon- assisted tunnelling provides 
a physically relevant model for which iV-th order time-dependent perturbation theory can be large 
although all previous orders are small. 
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FIG. 1. (Colour online) Sketch of l/4-'photon' resonance: 
One photon has enough energy to make four particles tun- 
nel (Huj = 4AE). The "photons" are time-periodic potential 
modulations in the kilo-Hertz regime. 



I. INTRODUCTION 

Optical lattices are an important system for research 
with ultra-cold atoms [TH3] . Experimental developments 
enable the creation of lattices of controllable double- well 
potentials [4] [5] that can be engineered such that the 
tunnelling between neighbouring double wells can be dis- 
carded [6], allowing treatment of the system as a single 
double well. Loading the lattice from a Mott-insulator 
state allows deterministic population of fewer than six 
atoms in each well [TJ [8] . Combined with the ability to 
count the atoms in each well, this makes the double well 
system ideal for investigating fractional photon-assisted 
tunnelling via periodic shaking of the lattice that is typ- 
ically a small effect in other systems [9UTTj. The "pho- 
tons" are time-dependent potential modulations in the 
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kilo-Hertz regime; a sketch of the 1/4-photon resonance 
can be seen in Fig. [I] 

Research on periodic shaking has focused on ef- 
fects ranging from destruction of tunnelling [12-14J over 
tunnelling-control [T5HT7] , field-induced barrier trans- 
parency [18], two-dimensional solitons [19], super Bloch 
oscillations [20-22 , phase-jumps (23J [24] and dynamics 
of bound pairs in optical lattices [25] [26] or NOON- 
states [27]. Complementary studies of transport of 
ultra-cold atoms in lattices include controlled trans- 
port of Bose-Einstein condensates (BECs) between two 
wells [28J. Experiments include control of superexchange 
interactions [29, 30 and directed transport via a Hamil- 
tonian quantum ratchet [3T] . 

While integer photon-assisted tunnelling [32j [33] is 
essentially a single particle effect that survives inter- 
actions, fractional photon-assisted tunnelling is a true 
many-particle effect which only occurs for interacting 
particles. Rather than being the small effect predicted in 
Refs. [SHU], fractional photon-assisted tunnelling like the 
1/2-photon resonance can be a large effect for two par- 
ticles per double well [30j [34] . A related experiment on 
photon-assisted tunnelling of strongly correlated atoms 
can be found in Ref. [35] . 

In this Letter we investigate fractional-photon assisted 
tunnelling of ultra-cold atoms in a double well poten- 
tial. The paper is organised as follows: Section [IT] in- 
troduces the model describing atoms in the periodically 
shaken lattice, whilst Sec. |III| explores the effect for both 
few and hundreds of atoms per well. In Sec. |IV| we 
demonstrate that for four particles, both the 1/4- and the 
1/2-resonance will provide clear experimental signatures. 
Section |V| explains a feature in the photon-assisted tun- 
nelling plot reminiscent of avoided crossings. In Sec. VI 
we show that for some parameters, even though the first 
few orders of time-dependent perturbation theory might 
be small, higher order perturbation theory can still cor- 



2 



rectly predict the fact that fractional photon-assisted 
tunnelling is a large effect. 



II. MODEL 

A. Hamiltonian 

For both the case of a few ultra-cold atoms, or a small 
BEC loaded into the double-well potential modulated at 
frequency lj, the system can be described using a many- 
body Hamiltonian with a two-mode approximation as fol- 
lows 033 Ez] 

H = (c\c 2 + 4&L ) + UK, (cjcfoc! + c\c\c 2 C^j 

+ fi(/i + Mi sm(cji)) (clc 2 - c\c^j . (1) 

The operators c -/ct annihilate/create a boson in well j; 

is the tunnelling splitting, H/jlq denotes the tilt be- 
tween well 1 and well 2 and h\±\ is the driving ampli- 
tude (cf. Fig.0. The on-site pair interaction is denoted 
by 2Hk. For calculations beyond this model see, e.g., 
Refs. [38H4T] . Here, we use the two- mode approximation 
as the experiment [30 demonstrates that this approxima- 
tion describes the physics of the 1/2-photon resonance. 
Our focus lies in identifying and understanding interest- 
ing signatures of photon-assisted tunnelling rather than 
quantitative predictions calculations including effects of 
higher energy levels (cf. Ref. [30 ) will subsequently de- 
pend on precise experimental details like the depth of the 
lattice. 

In order to characterise the photon-assisted tunnelling, 
we use the experimentally measurable time-averaged par- 
ticle transfer probability, 



which is based on the interaction picture, turned out to 
be useful to solve the Schrodinger equation [TO] , The 
Schrodinger equation then is equivalent to the set of N+l 
differential equations: 

iha v (t) = (z/li^+^M^K+iW (5) 
+ v = 0, 1, . . . N 

which uses the notation a_i(t) = a/v+i(£) = 0; the phase 
factors read: 

K(t) = (6) 
exp (i [2(N - 1 - 2v)Kt - 2fi t + 2/ii cos(ut) /u}) . 



B. Time-dependent perturbation theory 

In the following, we always use the experimentally re- 
alistic initial condition that all the atoms are in the lower 
well at t = 0, i.e. ao(0) = 1 [3 E]. Then at a later time 
t, zeroth-order time-dependent perturbation theory gives 
the Fock-state amplitudes as: 



4 0) (*) = i 

af\t)=0, j>0. 



(7) 

(8) 



The first non-zero order of is obtained for k = j: 



W(t)si #a 

J 2 



xJdTh j - 1 {Tya'fj: i 1) {T), j>l, (9) 



where * denotes the complex conjugate. 



C. Tunnelling dynamics 
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V(t))dt. (2) 



Rather than trying to understand the physics by di- 
rectly solving the Schrodinger equation corresponding 
to the Hamiltonian (IT), we derive an equivalent set of 
differential equations 34 . We choose the Fock basis 
\v) = \N — z/, v) for which the label v = ... N refers 
to a state with TV — v particles in well 1, and v parti- 
cles in well 2. In this basis, the Hamiltonian ([T]) can be 
written as the sum of two (N + 1) x (N + l)-matrices, 



H = Ho(t)+H! 



(3) 



where the diagonal matrix Ho(t) includes both the in- 
teraction between the particles and the applied potential 
difference while the non-diagonal matrix Hi contains the 
tunnelling-terms of Eq. ([T]) . 
The ansatz 



■ J (u\H (t)\u)dt 



(4) 



To understand the tunnelling dynamics, we expand the 
oscillatory term at frequency u in the phase factors of 
Eq. (I9|) in terms of Bessel functions [42] 



a iz cos(u)t) 



J2 Mz)i k e 1 



kujt 



(10) 



Including all these terms in analytic calculations does, in 
principle, lead to analytic results for the tunnelling. How- 
ever, evaluating these analytic formulae is numerically 
much more intensive than solving the time-dependent 
Schrodinger equation [10 . Combining the rotating- wave- 
approximation based approach which includes only the 
slowly oscillating terms of hj [34] with the above time- 
dependent perturbation theory leads to a simpler form 
for the non-zero perturbations of Eq. ([9| : 
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FIG. 2. (Colour online) Two-dimensional projection of the 
time-averaged particle transfer ^ as a function of both driv- 
ing frequency uj and interaction ^ for iV = 100 particles. 
The averaging time is QT — 100, the driving amplitude 
2/ii/cj = 1.8 and the tilt fio = The one-photon reso- 

nance is clearly visible: it starts at the point defined by k = 
and uj — and moves towards lower frequencies for larger 
interactions. While many-integer-photon-resonances like the 
two-photon resonance near k = and uj = 1.5Q are also 
visible, fractional photon-assisted tunnelling like the one-half 
resonance near uj — 6Q (cf. Ref. 9 ) are not visible on this 
scale because they are only a small effect for large particle 
numbers. 



where the integer kj is chosen such that \rjk\, with 

ril j) = -ku + 2/i - 2[N - (2j - 1)]« , (12) 

j = l...N. 

is minimised at k = kj. As in Ref. [34], it might some- 
times be preferable to minimise | Yjl I rather than 

each rj^ separately. For other cases (e.g. near a zero of 
one of the Bessel functions) more than one term will have 
to be included in the above sums. We define the number 
of "photons" involved in the tunnelling process, 



^photons 



N 

£ 



k n , 



(13) 



such that it corresponds to the total "energy" transferred 
(rather than taking the sum of the moduli). 



III. BEC OR FEW PARTICLES PER DOUBLE 
WELL? 

In order to demonstrate why it is necessary to use a few 
atoms per double well as opposed to a BEC to observe 
fractional resonances, we compare the results obtained 
for N = 100 particles with N = 4. Figure [2] shows the 
time averaged particle transfer probability for N = 100 
atoms initially in the lower well. While there are no 



FIG. 3. (Colour online) Time- averaged transfer probability 
to well 2 for N = 4 calculated using the same parameters as 
in Fig. [2] but with an averaging time of QT = 10. Experi- 
mentally, this could be realised in the double-well lattice of 

Refs. ram. 



resonances visible at higher frequencies, the one-photon- 
resonance which starts at uj = 3Q for small interactions 
is clearly visible. For a BEC initially in the upper well, 
the one-photon resonance would move toward higher fre- 
quencies for increasing interaction. 

Figure [3] shows the time-averaged transfer probability 
for N=4 calculated using the same parameters as Fig. [2] 
with an averaging time of QT = 10. For this shorter 
averaging time-scale, Fig. [3] displays many features also 
seen in Fig. [2] The visible lines are either tunnelling res- 
onances which can be understood on the single-particle 
level like the 1-photon resonance near k « and uj ~ 317 
or the horizontal lines. These correspond to, e.g., the 
energy of all particles being in the lower well and one 
particle having tunnelled being equal. The straight lines 
with non-zero gradient correspond to adding one or sev- 
eral photons to those horizontal lines. 

More interesting features, including fractional photon 
resonances, emerge for larger averaging times (Fig. [4|. 
While the short-time effects visible in Fig. [3] can be ex- 
plained by simply looking at tunnelling of a single parti- 
cle, Fig. [4] displays many features which are due to many- 
particle tunnelling. For larger particle numbers, similar 
features will only be visible on even larger time-scales 
(cf. Sec. |Vl|. In Fig. |I| both the 1/2-photon and the 
1 / 4-photon resonance are visible for a broad range of in- 
teraction strengths. 

Integer-photon resonances essentially are single- 
particle effects which survive interactions. Fractional 
photon assisted tunnelling, however, is a true many- 
particle quantum effect. For four particles, one ex- 
pects [11 to observe the 1/2- and the 1 /4-resonance both 
of which are clearly visible in Fig. [4] Odd fractions 
like the 1/3-photon resonance (cf. Sec. fvT| ) only occur 
for odd particle numbers N > 3 (cf. Ref. [H]). Frac- 
tional resonances also appear for the case of N = 100 
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FIG. 4. (Colour online) Time- averaged transfer calculated for 
the same parameters as Fig. [3] except for OT = 100. While 
the 1/2-photon resonance is clearly visible near uj — 60, the 
1 /4-resonance can also be found for some interactions near 
uj 120. 



(cf. Ref. [9]). However, for the experimentally motivated 
comparatively short timescales used both here and in 
Ref. [9], fractional photon-assited tunnelling is a very 
small effect even for small BECs. We therefore focus 
on the case of few-atoms per well for the remainder of 
the paper. 



IV. THE 1/2- AND 1/4- RESONANCE 

For fractional resonances in a double well lattice with 
few atoms per double well, most of the physics of the 
tunnelling process can be understood by taking the ap- 
proach of Sec. II C to a level beyond perturbation the- 
ory. Within the approximation motivated by the ro- 
tating wave- approximation, this leads to time-dependent 
(iV + l) x (iV + 1) matrices which can, in some cases, even 
be solved analytically (see the appendix). Figure [5 shows 
that the 1/2-photon resonance will be clearly visible for 
four particles. As the full width at half maximum is an 
order of magnitude larger than 1% (the typical error [43] 
with which the interaction can be fixed in experiments 
as Ref. [8]), the 1/2-photon resonance could be observed 
with the existing experimental setup of Refs. [SJ [30] . 

Given the fact that the 1/2-photon resonance has al- 
ready been observed experimentally [30] , it would be even 
more interesting to investigate the 1 / 4-photon resonance 
as N = 4 is the lowest number of particles for which it can 
be observed. Just because 4 particles produce a large ef- 
fect in the numerics does, however, not automatically im- 
ply that it is a true 4-particle effect: Many features which 
are already visible for the short averaging times of Fig. [3] 
would also occur at least similarly for lower particles. 
Thus, even if we restrict our search for a 1 /4-resonance 
to parameters near uo = 120, this could coincide with 
large tunnelling for lower particle numbers. Due to the 
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FIG. 5. (Colour online) 1/2-photon resonance for four par- 
ticles; displayed is the transfer to the second well at time 
t* = 13.115 as a function of the scaled interaction parame- 
ter k/Q. Dotted/magenta: model discussed in the appendix 



[Eq. (A14)], dashed/blue: full numerics. The parameters are: 
u = 60, fj,o = 1.50, 2/ii/O = 4.567. As the full width at 
half maximum is an order of magnitude larger than the typ- 
ical experimental accuracy for k/Q [43], the 1/2-photon res- 
onance could be observed with existing experimental setups. 
Choosing the correct time will also be feasible: For k — 1.50, 
numerics shows that missing the time t* by 10% only leads 
to deviations of Ptrans(t) from Ptrans(t*) ~ 0.92 by less than 
3%. 
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FIG. 6. (Colour online) a) 1/4-photon resonance for four par- 
ticles (« = 2.250, u = 120; 2/xi/cj = 3.08, /i = 1.50). Black 
line: exact numerics. The simplified model (magenta/grey 
line) correctly describes the dynamics for not too large times. 
In order to demonstrate that this is really an effect which 
only occurs for four particles, N < 4 was also investigated: 
b) N=3, c) N=2, d) N=l. Note that the scale on the vertical 
axis differs by an order of magnitude between a) and b)-d). 



experimental way to load double- well lattices via a Mott- 
insulator [30 , the harmonic confinement will prevent ex- 
periments with all double-wells being filled with exactly 
4 particles. However, as can be clearly seen from Fig. [6] 
only the wells initially loaded with 4 atoms can contribute 
to the observable signature of the 1/4-photon resonance. 
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FIG. 7. (Colour online) Two-dimensional projection of the 
fourth root of the time- averaged transfer for N — 2 particles. 
The averaging time is Tft = 500, the tilt is /io = 1.5ft and the 
shaking amplitude 2/ii/cj = 5.5201 (for which | Jo (2/ii/u;) | < 
IO -5 ). In Sec. [v[ these parameters allow the explanation of 
an example how the apparent avoided-crossing-type features 
(cf. Fig. [4| can appear. 
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Combined with the fact that r]\ 
line 
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3ft 



corresponds to the 
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this explains why transfer becomes large when following 
this line away from the point (uj = 6ft, k = 1.5ft); for 
the line perpendicular to this line through the same point 
such increased transfer is neither to be expected nor does 
it show in the numerics of Fig. [7] 



VI. 1/iV-RESONANCE: LARGE EFFECT ONLY 
IN iVTH ORDER PERTURBATION THEORY 



V. AVOIDED-CROSSING-TYPE FEATURES 

Figure [4] shows several features which resemble avoided 
crossings. The smallest particle number for which they 
can occur is N = 2 for which the feature is particu- 
larly strong at the 1/2-photon frequency near zeros of 
the Bessel function responsible for the tunnelling of the 
first particle. Such a situation is depicted in Fig. [7[ For 
parameters near the centre of this figure, i.e. /io = 1.5ft, 
uj = 6ft and k = 1.5ft, the tunnelling of the first parti- 
cle would normally be described by a 0-photon process 
while the tunnelling of the second particle would be a 
1-photon process. Only the average number of photons 
per tunnelling process justifies the word "1/2-photon res- 
onance". If, however, Jq (2/ii/a;) = (as chosen for 
Fig. [7]), we have an entirely different situation. Now 
the tunnelling of the first particle consists of two com- 
peting processes: a 1-photon process and a —1-photon 
process, making the overall tunnelling a superposition of 
a 0-photon and a 2-photon process [see Eq. (13)]. As 
J-i(x) = — Ji(x), we have: 



ft 

"71 
ft 

"71 



exp 



(irj^t^j exp (irj^lt^J 



(1) 



(1) 

iry_[ 



■ (1) ■ (i) 
irj\ ir]_{ 



(14) 



This equation already explains why there is no transfer 
visible at /io = 1.5ft, uo = 6ft and k = 1.5a;, as 



(15) 



for these parameters. In the vicinity of this point, the 
time-independent part of Eq. (14) will be responsible for 



The system investigated here offers the unique possi- 
bility to construct physically relevant examples for which 
the perturbation theory is small up to (TV — l)th order 
while the TVth order produces results that dominate the 
dynamics. 

If 



for all j = 1 ... N and furthermore [44] 



E 



= 



£<j<N 
1=1, j = N 



(18) 



(19) 



perturbation theory will give only (small) oscillatory 
terms up to (TV — l)th order. However, in TVth order 
we obtain a term linear in time proportional to: 



N 



leading 



n ^ ( 



(20) 



However, just because Eq. (20) becomes large this does 



ci2 eventually becoming larger, in particular, for param- 



not automatically imply that this can be observed in 
the numerics: as soon as a^\t) becomes large, ajy-i(t) 
changes which in turn influences ft at (t) • Nevertheless, 
there are examples for which Eq. ( |2Q| ) well describes the 
tunnelling: 

For N = 3, the above criteria are fulfilled, e.g., for 
/io = 1.5ft, n = 2.1ft, uj = 9ft and 2/ii/cj = 0.75: 



(21) 
(22) 
(23) 



The tunnelling can thus be understood to be a 1-photon- 
process, the photon being responsible for the tunnelling 
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FIG. 8. (Colour online) Probability to find n particles in 
the upper well if initially all N — 3 particles are in the lower 
well for jio — 1.50, k = 2.1Q, uj — 9Q and 2/ii/cj = 0.75. a) 
Black solid line: numerics using the model ([!]); magenta/grey 
dotted line: as predicted by numerically evaluated perturba- 
tion theory by adding |a^| 2 + la^l 2 as defined in Eq. 
The sum is much smaller than the probability to transfer all 3 
particles which is displayed in the lower panel, b) Solid black 
line: numerics using the model dotted magenta/grey line: 
la^l 2 as defined in Eq. ([9]), dashed blue/black line: analytic 
estimate of Eq. (24). 
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FIG. 9. Probability to find n particles in the upper well 
if initially all N = 4 particles are in the lower well for 
fio — 1.50, uj — 12Q k — 1.950, and 2/xi/o; — 0.5. a) Black 
line: numerics using the model |l]); magenta/grey line: as 
predicted by numerically evaluated perturbation theory by 
adding |a^ 1} | 2 + | 2 + |4 3) | 2 as defined in Eq. The sum 
is much smaller than the probability to transfer all 3 particles 
which is displayed in the lower panel, b) From bottom to top: 
numerics using the model |l]), analytic estimate of Eq. (29) , 



numeric evaluation of la^ 4 ^ 2 as defined in Eq. (|9|. The upper 



two curves lie very close together. 



of the third particle with no photons being involved in 
the tunnelling of the other two. Thus: 



(4 3 V)) 



3sr 



leading 4^> (t,™+t,™) 



(24) 



Figure [8] shows that Eq. ( 24 ) well describes the time- 
scale on which the probability for all particle having tun- 
nelled to the other well becomes large. The figure fur- 
thermore shows that the probability to find either one or 
two particles in the upper well remains small and this ef- 
fect can thus be labelled co-tunnelling of three particles. 
Equation (24) does, however, not precisely match the re- 
sult of the numerical perturbation theory of Sec. |IIB| 
Higher order processes would change the prefactors in 
both equations. For N = 3 there is, e.g., the +1-, — 1-, 
-fl-photon process which also contributes to the linear 
time-dependence, although it is low compared with the 
dominating 0-, 0-, -fl-photon process included in the an- 
alytic result. Note that changing the driving amplitude 
might change which is the dominating contribution. 

For N = 4 similar parameters can be found. Using 
/iQ = 1.50, uj = 120, n = 1.950, and 2\i\juj = 0.5, one 





= -8.70 




= -0.9O 




= 6.90 




= 2.70, 



with 



has as the leading order contribution: 

(25) 
(26) 
(27) 
(28) 

t . 

(29) 

This analytic function correctly predicts that the tun- 
nelling of all 4 particles at once takes place at a much 
longer time-scale than for three particles (Fig. [9]). As 
for the previous figure, the tunnelling can be labelled a 
co-tunnelling process, now of four particles. 
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VII. CONCLUSION 

Both perturbation-theory and beyond-perturbation- 
theory approaches can be applied to understand the 
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physics of fractional photon-assisted tunnelling in optical 
super-lattices. In some cases even analytical calculations 
are possible. While the 1/2-resonance has already been 
observed experimentally, we predict that both the 1/3- 
photon resonance and the 1/4-photon resonance should 
be observable experimentally. 

For the 1/3-resonance, an interesting phenomenon oc- 
curs: although both first and second order perturbation 
theory are small, third order perturbation theory never- 
theless correctly describes it is a large effect. This effect 
can also be found for larger particle numbers. It will, 
however, be much more realistic to try and observe this 
effect for 3 particles than for much larger particle num- 
bers: Both the time-scales for the transfer to the upper 
well becomes large and the model (fil) becomes less valid. 



In the same way in which it is possible to describe the 
1/2-resonance for some parameters alternatively as co- 
tunnelling of a pair of atoms [30 , we thus numerically 
observe co-tunnelling for both 3 and 4 atoms. 
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for | a at | to be large 



7\. Furthermore, | being large is, 
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Appendix A: Understanding the tunnelling for N — 4 
via beyond-perturbation theory calculations 

Reference [34 demonstrates that it is possible, at least 
for cases with one dominant frequency for each tunnelling 
event, to use a beyond-perturbation theory version of 
Sec. |II C[ In the present Section, the analytic calculations 
within the effective model [34] are performed in detail for 
N = 4 particles. 

Using the abbreviations 



y/jy/N-j + l 



B 



k ' 



nj k (2^/uj) (Al) 
(A2) 



we can rewrite the equations of Sec. |IIC| in such a way 
that they can be solved beyond the perturbation theory 
approach. We are, however, restricted to cases where 
only one frequency plays a role. For those we have: 



d_ 
dt 



( ao(t) \ 
ai(t) 
a 2 (t) 
a 3 (t) 

\ a 4 (t) } 



= A(t) 



/ «o(*) \ 

a 2 (t) 
a 3 (t) 
V o 3 (t) / 



(A3) 



where A(t) is a Hermit ian 5x5 matrix which has only 
non-zero entries for (A(t)) n>m with n = m ± 1: 

(A(i))„, n+1 = B ( k n J exp (-ivjtft) (A4) 



and 



(A(t)) n+1 , n = (B^)*exp(4^). 
This can be solved using the ansatz 
/ 



(A5) 



fa (t)\ 
ai(t) 
a 2 (t) 

\ a 4 (t) J 



ao exp(— iujt) 

ai exp[-i(u-rj^)t] 
r •( C 1 ) (2)wi 



n {1) 

~ r v (!) ( 2 ) ( 3 )wi 

a 3 exp[-z(cj - ri K k ' - 77^ - r) K k J )t] 

\ a 4 exp[-z(cj - rj ki - f] k2 - f] ks - rj k4 )t\ / 

(A6) 



which leaves us to calculate the eigenvalues of the time- 
independent matrix: 



/ 



V 















-< 





d*(2) 
n k 2 


(1) , (2) 








d*(3) 

^3 















-5 



(3) 



CO 



«(4)- 

/C4 







-5 



(4) 



-5 



/C4 

(A7) 



In order to obtain a 1/2-photon resonance the following 
conditions [10] are imposed on the choice of the ideal 

frequencies rj k ^: 



(2) 
77 (1) 



(4) 
(3) 

(2) 
% 2 



&2 



&4 



1, 
1, 



-^(fci + fc 2 )w-^o- (A8) 



Furthermore cj/2 = 2/io and thus one may choose ft = 
uo/A implying k\ + k 2 — — 1. Once choosing fci all fre- 
quencies are fixed. With the above conditions fulfilled 
it is now possible to solve the eigenvalue problem ana- 
lytically. For the sake of readability from now on the 
following further changes in the notation will be adopted 
in the four-particle case: 



n£> -> fix 



"fei+2 



(4) 
l-fel 



^3 
fi 4 



(A9) 



The eigenvalues are: 



^1 = 0, 
^2/3 = ^ I ^1 T + 2 



^4/5 = o 1 Vi T +2 



(Of + fl 2 + ^2 + ^2 _ I 

(n? + o 2 , + n§ + ol + y/itj I 



(A10) 



where K represents the expression: 
K = Of + Sl\ + 20 2 ,(0 2 - 



2n?(n| - n§ - n2) + (fi§ + ^I) 2 • (ah) 



The cor resp onding eigenvectors are inserted into the 
ansatz (A6) and yield five independent, not yet nor- 
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malised, solutions to the time-dependent problem: 



(*) 



,(3) 



(*) 



^4 


\ 






















1 


J 



2J7o ^3 ^4 



2172^3^4 



-k'+i 



-Vk) 



n 4 



ir)\t 



1 



-iu>2t 



Qi(-^i-^2+^3+^4 + ^^) 



2^2^3^4 



2^2^3^4 
i(^i+^2+^3 



2^73^4 

-fc'+l 



(*) 



(*) 



2^2 ^3 ^4 



2^2 ^3 ^4 



21^3174 
04-^e 
1 



.£771 1 



21^2 ^3 ^4 



21^2 ^3^4 



-Ql + y/K) 



2^3 r^4 



/ 



nl±v^) W5eiJ)1 t 



(A12) 

Any state of the four-particle system can now be 

expressed as a linear combination \^{t)) = biv^(t) + 
b 2 v^\t) + b 3 v( 3 \t) + M (4) + M (5) of the eigenstates. 
The coefficients 6^ are again obtained from the initial 
condition |tf(0)) = |0). 

The probability to find the system in a Fock state \i) 
is given by the overlap with \^(t)) and may be expressed 

as the square of an amplitude |<^(t)| 2 with being the 
i-th component of the j-th eigenvector: 

2 



\^{ t )\i)\ = \ ai {t)\ = 



i=i 



v i+i 



(*) 



\vU)(t)\ 



(A13) 



Finally, the time-dependent transfer Ptrans(^) is expressed 
in terms of the amplitudes a^(t) as 



4 |a 4 (i)| 2 + 3 |a 3 (i)| 2 + 2 |a 2 (i)| 2 + K(t)| 2 



(A14) 

This expression for given time and amplitude of the driv- 
ing is displayed as a function of the coupling strength in 

Fig m 



